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Abstract 

We use gauge/gravity duality to explore strongly coupled superconductors with dynamical exponent z = 
2. In the probe limit we numerically establish background solutions for the matter fields and plot the 
condensate versus the dimensionlcss temperature. We then investigate electromagnetic perturbations in 
order to compute the AC conductivity and also calculate the spectral function. Our results for the condensate 
, and conductivity are qualitatively similar to those of the AdS superconductor. However, we find that 



(for both s- and p-wave) the condensate does not approach a constant at very low temperature and the 
conductivity goes to one from below but never exceeds it in the high frequency limit, in contrast to the 
AdS black hole. We do not see a peak at nonzero frequency in the imaginary part of the AC conductivity 
along the x direction for the p-wave case. These features are due to the nontrivial dynamical exponent. 
To be specific, the black hole geometry considered in this work is anisotropic between space and time, very 
different from the Schwarzschild-AdS black hole, which results in different asymptotic behaviors of temporal 
and spatial components of gauge fields than those in the Schwarzschild-AdS black hole. 
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I. INTRODUCTION 



The application of gauge/gravity duality j| to the investigations of strongly coupled systems 
has gained broad interest ranging from QCD phenomena at low energy to strongly correlated con- 
densed matter physics, see, e.g., [§] for some recent reviews. Holographic superconductors have 
been constructed in || [|] by putting the Abelian Higgs model or SU(2) gauge field into the AdS 
black hole spacetime. When the Hawking temperature is decreased to a critical value, the black hole 
becomes unstable against small perturbations and develops hair by condensing some field to stabi- 
lize the system. This can be considered as the holographic realization of the superconducting phase 
transition. This kind of construction for holographic superconductors takes the (asymptotically) 
AdS black hole spacetime as the starting point. According to the AdS/CFT correspondence, the 
AdS black hole geometry corresponds to a relativistic CFT at finite temperature. However, many 
condensed matter systems do not have relativistic symmetry and it is therefore very natural and 
interesting to generalize these holographic superconducting models to non-relativistic situations. 

On the other hand, inspired by the dynamical exponent in condensed matter physics near 
the critical point, many papers have appeared on the construction of black hole geometries with 
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anisotropic scaling, such as the Lifshitz black hole. The dual geometry of Lifshitz fixed points was 
first proposed in || and then generalized to finite temperature in Q. In short, we now have a 
geometrical realization of a strongly coupled anisotropic field theory at finite temperature. It is 
expected that non-relativistic AdS/CFT will help us to understand some puzzles in unconventional 
condensed matter physicsQ. 

In H a variety of strange metallic behaviors have been realized by using gravity duals of the 
Lifshitz fixed points and some string theoretical realizations of this geometry have been proposed in 
that paper^]. More recently, holographic fermions have been studied in |l0| to produce a non-Fermi 
liquid behavior. With these studies in mind, we now generalize holographic superconductors under 
the relativistic AdS/CFT framework to the Lifshitz black hole geometry in order to explore the 
effects of the dynamical exponent and also in the hope to distinguish some universal properties of 
holographic superconductors. In actual fact, holographic s-wave superconductors with a Lifshitz 
fixed point have been constructed in [|ll], 12] and in Hofava-Lifshitz gravity in 113]. However, the 



work of [12 1 produced only the condensate and the results of [11| seemed surprising, especially the 
very small real part of the AC conductivity. We therefore more carefully study these systems here 
in order to clarify some confusions which appeared in [11] and also reveal some properties of p-wave 
superconductors with Lifshitz scaling. 

Although there is a dynamical exponent, which makes the Lifshitz geometry behave quite dif- 
ferently from asymptotically AdS spacetime, we find that qualitative behaviors of holographic 
superconductors with Lifshitz scaling (for both s-wave and p-wave) are basically the same as those 
of the AdS black hole case. The condensate has mean-field behavior near the critical temperature. 
A gap does appear when decreasing the temperature, which is manifest once looking at the real 
part of the AC conductivity. There is also a delta peak near zero frequency for the conductivity, 
which is a signal of DC superconductivity. The real part of the conductivity approaches one in the 
high frequency limit where the imaginary part goes to zero. These common characteristics of holo- 
graphic superconductors appear to be robust phenomena and can be taken as universal properties 
of gauge/gravity duality when applied to the study of strongly coupled condensed matter systems. 

We also find some other interesting features in our work. The first is that all the condensates do 
not approach some constant in the zero temperature limit compared to the BCS superconductor 
and the AdS black hole holographic superconductors. The real part of the conductivity never 
exceeds one and the imaginary part of the conductivity just approaches zero from above but never 

1 Actually, there appear many works by taking this kind of non-relativistic metrics to reveal some strange features 

of condensed matter system, see |Q] for an incomplete list. 

2 Another string theory realization of Lifshitz-like fixed points was investigated in B . 



goes below zero in the high frequency limit. What is more striking is that there is no pole at 
nonzero frequency for the imaginary part of the AC conductivity cr xx (u)) in the p-wave case. We 
therefore attribute these nontrivial features to the effect of the Lifshitz scaling. More specifically, 
the black hole geometry considered in this work is anisotropic between space and time, very different 
from the AdS black hole, which results in different asymptotic behaviors of temporal and spatial 
components of gauge fields than previous conclusions in AdS black hole. 

Our paper is organized as follows. In section [ll| we briefly review basic aspects of asymptotically 
Lifshitz black holes for further studies. Section |l| is concerned with the construction of s-wave 
superconductors in the Lifshitz black hole geometry. With the numerical solutions one can plot the 
condensate as well as the free energy difference between the normal and superconducting phases 
versus the dimensionless temperature. We find that the free energy difference is always greater than 
zero when below the critical temperature T c , which proves that the superconducting phase is at 
least thermodynamically stable. Then we move on to investigate the electromagnetic fluctuations 
of the system and numerically calculate the AC conductivity using linear response theory. The 



spectral functions for the electromagnetic perturbations are calculated as well. In section |V| the 
corresponding results for the p-wave superconductor are presented. A short summary is given in 
section |V[ 

II. HOLOGRAPHIC SETUP: THE GRAVITY DUAL OF THE LIFSHTIZ FIXED POINT 

In this section, we concisely provide some backgrounds for the gravity dual of the Lifshitz 
fixed point. As mentioned in section |, there exist field theories with anisotropic scaling symmetry 
between the temporal and spatial coordinates. This is found for example in some condensed matter 
systems near the critical point, 

t X z t, x i Xx\ (1) 

where z is called the dynamical exponent. One geometrical realization of this scaling symmetry 
comes from the generalized gauge/gravity duality: we can map this scaling symmetry in the field 
theoretical side to some geometrical symmetry in the gravity side. Then it is straightforward to 
write down the metric with this type of scaling symmetry, 

ds 2 = L 2 (^-r 2z dt 2 +r 2 J2 dx 2 + ^ j , (2) 

where < r < oo and L is the radius of curvature of the geometry. This geometry was first 
proposed in where the action sourcing this geometry was also given. The scale transformation 



takes the following form, 

t -> X z t, x i -> Ax*, r^ T -. (3) 

A 

When z = 1, the above geometry reduces to the usual AdSd+2 spacetime and the symmetry group 
is enlarged to SO(d + 1,2). 

According to the gauge/gravity duality, to put the dual field theory at finite temperature one 
can study the metric with a black hole and the Hawking temperature of the black hole is identified 
as the temperature of the dual field theory. It is therefore of great interest to construct a finite 
temperature version of the above geometry, i.e. the so-called Lifshitz black hole. However, it is 
difficult to obtain analytic black hole solutions in Lifshitz spacetimes. There are some attempts to 
construct Lifshitz black holes[] and we here follow the work of || . It was found that action of the 
form 

S = 16?r ^ +2 J d d+2 x^~g (r-2A- \d^<P - \e^T^ v T^ , (4) 

where A is the cosmological constant, is a massless scalar and Tu V is an abelian gauge field 
strength, admits the following black hole geometry, 

* 2 = L 2 (-r 2 'f(r)d? + r 2 £ dx\ + J^j , (5) 



To support the above black hole geometry, one also needs to give the backgrounds for <f> and J 7 , 



it ■ 



p \<t> - r ~ 2d \2 _ _2d_ 



T rt = q r z+d -\ q 2 Q = 2L 2 (z - l)(z + d). (7) 

Evidently, choosing the dynamical exponent z to be 1 reduces the Lifshitz black hole to the 
Schwarzschild AdS black hole in d + 2-dimensions. The Hawking temperature and the entropy of 
the black hole are, 

z + d z L d V d d 

where denotes the volume of the d-dimensional space. 



See Jl4|, ^i) for an incomplete list of the constructions of black hole geometries with anisotropic scaling symmetry. 
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We do a coordinate transformation u = vq/v to map the holographic direction r into a finite 
interval [0,1] as we find that it is more convenient to use this coordinate system when carrying out 
numerical calculations. What's more, we will focus on a 4-dimensional bulk theory and choose the 
dynamical exponent z = 2. With these choices, the bulk geometry is reduced to, 

ds 2 = L*(- r ^J(u)de+ r l(dx* + dy*) + 4^), f(u) = l-u^. (9) 

In this new coordinate system, the horizon is located at u = 1 and u = denotes the conformal 
boundary where the dual field theory lives. 



III. HOLOGRAPHIC S-WAVE SUPERCONDUCTORS WITH LIFSHITZ SCALING 



This section and the next are the central parts of our work. In this section we focus on aspects 



of s-wave superconductor. In subsection [III A| we list the equations of motion for the background 



fields and solve them by a shooting method. Subsection MB is devoted to the studies of the 
electromagnetic perturbations of the system. 



A. Solution for the background fields 

In jn]], it was shown that a charged AdS black hole supports charged scalar hair if the charge 
is large enough. Later this idea was generalized to the neutral AdS black hole in Q, which is the 
first model of an s-wave holographic superconductor. More specifically, it is constructed from the 
Abelian Higgs model in the AdS black hole background. The action for this system is, 

S = J^x^fR + ^-^F^-ld^-iqA^-VimY (10) 

where we choose y(|^|) = m 2 |\I , | 2 for simplicity. This is the minimal Lagrangian of the gravi- 
tational dual which holographically describes a superconducting phase transition. As mentioned 
before, we work in the probe limit, i.e, the black hole geometry is fixed and feels no effect of 
the matter fields. Above the critical temperature, the black hole background is stable and the 
scalar field \t can be set to zero. This corresponds to the normal phase. When the temperature is 
decreased to the critical value, the black hole background becomes unstable against small pertur- 
bations and the scalar field wants to condense in order to stabilize the system. Once this happens, 
the black hole develops hair and the system goes through a superconducting phase transition. 
The scalar field VP holographically models the order parameter of a conventional superconductor 
and gets a nontrivial profile only in the superconducting phase. As the AdS/CFT correspondence 
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maps a strongly coupled field theory to a weakly coupled gravity system, the holographic method 
is expected to give a description of strongly coupled superconductors in contrast to conventional 
BCS theory. However, it is still far from being clear as to the pairing mechanism in holographic 
superconductors. Due to the conformal characteristic of the AdS space, the chemical potential is 
usually introduced to explicitly break the conformal invariance and to make the temperature scale 
meaningful. This is achieved by turning on the time component of the U(l) gauge field. 

Having fixed the black hole geometry, the equations of motion for and \l/ can be easily 
extracted from the Euler-Lagrange equations and are listed here for further studies in later sub- 
sections, 

L =<9 M (V^gF^) = iq[ty*(d u $> - iqA u ^) - + iqA v **)], (11) 



d„[V=g(d^ - iqA^)} = ^[V'fltfD-jL + iqA^d^ - iqA^)}. (12) 



2|* 

From the radial component of eq. (|ll]) one can show that the scalar field can be taken as real. 
Therefore, the ansatz for the backgrounds of the gauge and scalar fields are 

A = (p(u)dt, * = V(«). (13) 

Recalling the Lifshitz black hole background given in eq. (||), we obtain the following equations of 
motion, 



„ z-1 . 2ib 2 

»" + </>' - -2jr-,<P = o, 14 



f'(u) z + 1 



,, u 2z 2 (p 2 , m 2 L 2 , . . 

f 2 {u) u 2 f{u) 



f(u) u 

where we have for simplicity chosen the scalar potential V{^) as m 2 ^ 2 , giving mass to \E', and set 
the charge q 2 L 2 = 1. The prime in these equations denotes derivative with respect to u and this 
notation will be used in the following. 



Before continuing, some remarks are in order. Actually, in eqs. (14,15), we have rescaled the 
fields (p and i/j for convenience. When the dynamical exponent z is nontrivial, the rescaling is 
different between ip and (j), 

The superconducting phase transition is a spontaneous breaking of the electromagnetic symmetry, 
which should be reflected in the holographic description in some sense. Strictly speaking, the gauge 



symmetry in the bulk corresponds to a global symmetry on the boundary filed theory according to 
gauge/gravity duality and therefore we have no electromagnetic symmetry on the boundary. Then 
the holographic superconductor model should at most be thought of as holographic superfluidity. 
However, as discussed in || this model can produce many features of superconductors and we can 
ignore this subtlety. 

We now have a look at the asymptotic behaviors of the background fields eft and ip. Near the 
horizon, one must have </>(l) = for its norm to be finite and the scalar field should also be finite 
there. Near the conformal boundary, we have the following asymptotic behavior from the Probenius 



analysis of eqs.(|Tj) and (15) near the singularity u = 0, 

(j){u -> 0) ~ p + pu 2 ~ z (z^2) or <f)(u ->• 0) ~ p + p\ogu (z = 2), (17) 

and 

V>(u -> 0) ~ ^ u +\ u + + ^ V" , (18) 

where the scaling dimension v± of the scalar operator O dual to the bulk scalar ip is given by 
v± = 2+2=l= v / ( g + 2 ) + 4m ^ j^. - g gjqjjjgjt that the BF bound for the scalar mass in Lifshitz background 
is now changed to m > — 4 for 4-dimensional bulk. To obtain explicit behavior for these two 
fields near the conformal boundary, we need to specify the mass square of the scalar field and the 
dynamical exponent. We here choose z = 2 and m 2 = 0, —3 (which are both above the BF bound). 

For the exponent z = 2, the boundary behavior of the 4>(u) field has a logarithmic termF] due to 
degenerate indices of eq. (|l4|) near the conformal boundary u = 0. Now the physical result is also 
changed: we should identify the constant term p as the charge density as it is a normalizable mode 
with respect to the logarithmic term and the non-normalizable logarithmic term p as the chemical 
potential. Fortunately, this happens only for the time component of the gauge field and the final 
result for the conductivity does not change too much. We will come to see this when studying the 
AC conductivity along the x direction for the p-wave case. 

From the boundary behavior of the scalar field ip, we can directly read off the expectation value 
of the dual operator O. When m 2 = 0, the indices v = and 4, 

ip(u->0) ~^ (0) +V (4 V- (19) 

To make the superconducting phase transition a spontaneous breaking of symmetry, we should 
impose that 

V> (0) =0 and (0 4 }~V> (4) - (20) 

4 This kind of asymptotic behavior for the bulk field also happens in the study of p-wave holographic superconductor 
by taking the D3/D3 model in CI. 



Therefore, there is only one theory for this choice of the scalar mass squared. 
While for m 2 = —3 the indices v = 1 or 3, 

ip(u ->■ 0) ~ i/) (1 Wi/) (3) m 3 . 



(21) 



Evidently, the two modes in the above equation are both normalizable according to [18|, and in 
order to make the theory stable, we should either impose 



(i) 



and (0 3 ) ~ V 



(3) 



(22) 



or 



V> (3) = and (0i ) ~ ij> 



(i) 



(23) 



We now have two theories for m 2 = —3 corresponding to dimension 1 or dimension 3 order param- 
eters, respectively. In the following analysis, we will concentrate on the dimension 3 case because 
we find that the ground state of the dimension 1 theory is numerically unstable. 

With the boundary conditions mentioned above, we can now use a numerical shooting method 
to solve the coupled nonlinear eqs. (|TJJl^). The condensates corresponding to operators O3 and 
4 are plotted in FIG. |. 
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FIG. 1. The condensates of the s-wave superconductor for the dimensionless 3D operator O3 (blue): 



{/(0 3 )/T c and 4D operator O4 (red): {/(0 4 )/T c 



From this figure, we see that the condensates go to zero at the critical temperature T c . However, 
they do not approach some fixed constants as the temperature T — > 0, which is different from 
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both conventional BCS theory of weakly coupled superconductors and AdS black hole holographic 
superconductors. We may attribute this effect to the nontrivial dynamical exponent z 7^ 1. We 
also find that the expectation values for the operators y/ (O3) and \J (O4) are much larger than 
the BCS predictions at zero temperature. This is consistent with the results of AdS black hole 
holographic superconductors. Perhaps this is due to the fact that the strongly interacting nature 
of the holographic superconductor in contrast to the BCS theory. 

In the mean field theory for the superconductor, the order parameters have a square root 
behavior near the critical temperature T c , 

(O) ~ (T c - T) 1 / 2 . (24) 

By fitting these curves, we find the mean field behavior also holds in our results. Specifically, for 
the dimension 3 theory: 

(O3) « (18.9117T C ) 3 (1- T/T,) 1 / 2 when T -> T c , (25) 

where the critical temperature T c ~ 0.0351935/x. For the dimension 4 theory: 

(C 4 ) ~ (35.5940T C ) 4 (1 - T/T c ) 1/2 as T -> T c , (26) 

with the critical temperature T c ~ 0.0229931/U. 

The relation between the chemical potential and charge density is plotted in FIG. ^. We see 
that the profile has a small deviation from the linear behavior in the superconducting phase. There 
is one critical value for the chemical potential where the charge density becomes nonzero. This is 
actually the critical point where the superconducting phase transition occurs. 

At the end of this subsection, we plot the free energy for the two theories, which can be taken as 
one piece of evidence that the superconducting phase transition does happen when the temperature 
is decreased to the critical one. With the equations of motion for the backgrounds <f> and ip, we can 
reduce the action (??) to some simpler expression by integrating by parts. Then, the free energy 
difference between the normal and superconducting phases is 

A=„_ so _ ||,,.« £ iu^ )[m m? = ^<^F. (27) 

FIG. H is for the plot of the free energy difference. Clearly, the superconducting phase is thermally 
stable when T < T c . 
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10 20 30 40 50 60 70 

FIG. 2. The profile of the charge density as a function of the chemical potential in the superconducting 
phase for the s-wave superconductors: the blue one is for the dimension 3 theory and the red one is for the 
dimension 4 theory. 
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FIG. 3. The dimensionless free energy difference between the normal and the superconducting phases of the 
s-wave superconductor for O3 (blue) and O4 (red) theories, respectively. 

B. Fluctuation analysis: AC conductivity and the spectral function 

To calculate the AC conductivity of this system, we need to study the electromagnetic per- 
turbation of the above background (here, we do not consider fluctuation of the scalar field). For 
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the s-wave superconductor, the conductivity is isotropic and it is therefore equivalent to study 
any spatial component of the gauge field. We here choose radial gauge for the electromagnetic 
fluctuation, i.e., a u = 0. For the purpose of the AC conductivity, we take the following ansatz for 
the electromagnetic perturbation a x (t,x,u) 

a x (t,x,u) = e- iwt a x (u), (28) 

where we have taken spatial momentum to zero. Then, the equation of motion for a x (u) is, 



f{u) z - 1 
f(u) u 



a' + 



Co 2 u 2z ~ 2 2^ 



f{uf u 2 f(uf 



0, (29) 



where the dimensionless frequency Co is defined as Co = (z + 2)w/47rT. According to linear response 
theory, the conductivity is given by the Kubo formula, 

, , G(u,k = 0) , , 

ILO 

where the retarded Green's function G{u,k = 0) for the operator dual to gauge field can be 
computed according to the prescription given in |l9| . 

Near the horizon, we should take the ingoing wave boundary condition for the electromagnetic 
field fluctuation in order to calculate the retarded Green's function, 

a x {u) = (1 - n)-^/ 4 [l + al(l - u) + a 2 x {l - u) 2 + o^l - uf + • • • ], (31) 

where we have set the scale of a x to be one by making use of the linearity of the fluctuation equation. 
The coefficients in the above expansion can be uniquely determined once the background i/j and 
the frequency Co are specified. 

Near the conformal boundary u = 0, the general solution to the fluctuation is of the form, 

a x {u) = A° x + A l x u 2 + ■ ■ ■ . (32) 

The conductivity can be expressed as 

'<"> - hi < 33 > 

by using the Kubo formula (|30|). By the way, we also give the definition of the spectral function 
1Z plotted later, 

A 1 

K = -2ImG(w, k = 0) = -2-f . (34) 



Generically, the equation of motion for the electromagnetic fluctuation ( p9[ ) cannot be solved 
analytically due to the presence of the scalar field. The philosophy of the numerical method is that 

12 



we can use the power series solution as in eq. ( |3l"|) and do numerical integration from the horizon 
to the conformal boundary. Then we can extract the coefficients A® and A\ and therefore get the 
conductivity. Our numerical results for the AC conductivity are plotted in FIG. [| and FIG. |5|. 




FIG. 4. The real (imaginary) part of the AC conductivity of the s-wave superconductor for the dimension 
3 theory at T/T c = 1.0,0.798694,0.269798,0.0864351 from top to bottom (bottom to top). 




50 100 150 200 250 300 50 100 150 200 250 300 



FIG. 5. The real (imaginary) part of the AC conductivity of the s-wave superconductor for the dimension 
4 theory at T/T c = 1.0,0.565265,0.298191,0.0879648 from top to bottom (bottom to top). 

From the plots, we find that the real part of the conductivity approaches one in the high frequency 
limit, which is equal to the results of the normal phase. This is explicit when taking a look at 
eq. (p9|): the effect of the condensed field tp can be neglected at large Q. Compared to the result 



for the conductivity in [11], our result is of order 1 although these profiles are qualitatively similar. 
We found that the results for dimension 3 and dimension 4 theory are qualitatively the same. We 
also see that a gap forms as the temperature is lowered and the gap gets increasingly deep until 
the conductivity is exponentially small, which is the same as the AdS superconductor 

In the imaginary part of the conductivity, there is a pole at zero frequency. This can be explained 
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from the Kramers-Kronig relation 

ImW„)] = -P r ^5fM^, (35) 

where P denotes the principal value of the integration. Prom this formula, we can see that the 
real part of the conductivity contains a delta function, Re[cr(a;)] = it5(uj), only when the imaginary 
part has a pole, Im[<r(a;)] = 1/u. Actually, there is a peak at zero frequency for the imaginary part 
of the conductivity as can be seen in the right plots of FIGs. [| and [5||] 

In FIG. H we plot the results for the spectral function corresponding to the gauge field fluctuation 
a x . The two figures are consistent with the real part of the conductivity: a gap will appear when 
the condensate is nonzero (equivalently, when the temperature is low enough) and the spectral 
functions have a linear behavior with respect to the frequency at large u/T (which corresponds to 
large u/T behavior of the conductivity). 




T 

FIG. 6. The spectral functions of the s-wave superconductor for the dimension 3 theory 
(blue) at T/T c = 1.0,0.798694,0.269798,0.0864351 and the dimension 4 theory (red) at T/T c = 
1.0, 0.565265, 0.298191, 0.0879648 from top to bottom. 



In fact, we can explicitly plot the delta peak in the real part of the conductivity. However, when carrying out the 
numerical calculations, it is a little hard because we find that this peak will appear at about uj ~ 10 -20 and the 
numerical computations are not stable there. 
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IV. HOLOGRAPHIC P-WAVE SUPERCONDUCTORS WITH LIFSHITZ SCALING 



In this section we present the results of the p-wave superconductor in the Lifshitz black hole 
background. The structure of this section is the same as the previous one. 

A. Solution for the background fields 

As we have seen in the last section, the s-wave holographic superconductor is very simple and 
also elegant in describing some important features of superconductors. However, the Abelian- 
Higgs model for the s-wave superconductor appears to be less universal: we have to specify a 
potential for the scalar. On the other hand, the p-wave superconductor has already been observed 
in some experiments of condensed matter physics. It is natural as well as interesting to extend 
the construction of holographic s-wave superconductors to the p-wave situation. This has been 
achieved in [Q] by introducing the SU(2) gauge field into the AdS^ black hole background^. In this 
p-wave model the chemical potential and the order parameter have been unified into one field, the 
nonabelian gauge field, and the action is uniquely determined by gauge invariance. In this sense, 
the p-wave holographic superconductor is more universal than the Maxwell-Scalar system for the 
s-wave superconductor. The action is simply of the form, 

S = i / d ' x ^[R + - \f%F^\ (36) 

where the SU{2) gauge field strength is defined as F° v = d^A^ - d„A® + e^A^A^. The Yang-Mills 
equation for the gauge field is 

cy7^F a H + y/^e abc A b ^F c ^ u = 0. (37) 

The gauge field ansatz is taken as 

A = cf)(u)T 3 dt + ^(u^dx, (38) 

where r a is the generator of the SU(2) gauge group. Plugging this ansatz into the Yang-Mills 
equation ([37]) results in 

f + —4>' - = 0, (39) 

u f(u) 

6 Another approach to p-wave superconductor where the boundary field theory is known is based on D-brane probe 
in black brane geometry Eot. 
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and 



/(«) 



u 



,2z-2a2 



/(«) 



-V> = o. 



(40) 



As in the s-wave model, boundary behavior of the background fields <j>(u) and i/j(u) is listed as 
below having chosen z = 2, 



0(u — > 0) ~ p + /x log u, 
ip(u -»• 0) ~ V> (0) + ^ (2) " 2 - 



(41) 
(42) 



The logarithmic term appears again for the field, which suggests that we should identify the 
constant term p as the charge density. We also impose ^ (0) to be zero in our numerical calculations 
to insure that the superconducting phase transition is a spontaneous breaking of the symmetry. 
Then the value of is proportional to the expectation value of the p-wave order parameter, i.e. 
(O2) ~ ■ Our numerical results are plotted in the following figures. The free energy difference 
between the normal and superconducting phases is plotted in FIG. 0. As in the s-wave model, the 
superconducting phase is stable compared to the normal phase when T < T c . The charge density 
as a function of the chemical potential can be found in FIG. ||. Explicitly, these plots are very 
similar to the s-wave case. 
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FIG. 7. The free energy difference between the normal and the superconducting phases for the p-wave. 



The condensate of the p-wave order parameter (O2) is plotted in FIG. ||. It is clear to see from 
this plot that the condensate does not approach some fixed value when T — > 0, which is the same 
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FIG. 8. Density as a function of /j, for the p-wave superconductor. 

as the results of the s-wave model mentioned in the last section. Once again, this value is much 
larger than that of the AdS superconductor. By fitting the curve near the phase transition, we 
find a mean field behavior for the condensate, 

(0 2 ) ~ (17.8217T C ) 2 (1-T/T C ) 1/2 when T -> T c , (43) 

where the critical temperature is T c « 0.0367064/i. 
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FIG. 9. The condensate < C 2 > of the p-wave superconductor. 
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B. Fluctuation analysis: AC conductivities and spectral functions 



Owing to the appearance of the background for the gauge field, say Aq, the SU(2) gauge sym- 
metry is explicitly broken to the U(l)3 generated by the rotation in the colored 12 plane. This 
residual symmetry is identified as the electromagnetic symmetry. For the conductivity, we need 
to perturb the system and investigate its linear response. The conductivity of the p-wave super- 
conductor is anisotropic, which makes the fluctuation analysis more complicated than the s-wave 
case. We here focus on the decoupled mode a^(t, it), decoupled sector {a*(t, it), a 2 (t, it), a^.(t, it)} 
and leave the full fluctuation analysis including spatial momentum for future investigations. 

For the fluctuation mode a^(i,u) ~ e~ lu)t ay(u), which is decoupled from other modes, we have 
the following equation, 



4 + 



/fry 
/(«) 



4 + 



u 2z ~ 2 0? 



f(u) 2 f(u) 



(44) 



This mode looks like the electromagnetic fluctuation in the scalar system for the s-wave model and 
the procedure for numerical computation of the conductivity along the y direction is of course in 
parallel with the s-wave situation. We here only present the final results for the conductivity in 



FIG. 1C 
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Im(tr„) 20 




FIG. 10. The real (imaginary) part of the AC conductivity along the y direction of the p-wave superconductor 
at different temperature T/T c = 1.0,0.588587,0.171829,0.079889 from top to bottom (bottom to top). 



We then move on to the conductivity along the x direction. For this purpose, one should analyze 
the fluctuations {a|(t, u), a 2 (t, u), a^.(t, u)} ~ e~ lujt {aj (u), a 2 (it), a^(it)}, which is also decoupled 
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from other modes but is self-coupled. There are three coupled equations of second order, 



al + 



z - 1 



-ai + 



<N> „3 

/(«) 



o. 



3" 

ax + 



u /(„) 
/(«) u 



(iwa^ + i/>af) = 0, 



(45) 



,2z-2 



/(«) 



and two constraints of first orderQ, 



itDal' + <M 2 ' — <^>'a 2 = 0, 



•~ 2' / 1' 

iooa t + (pa t 



u 



(46) 



0. 



Near the horizon u = 1, we choose the ingoing wave boundary condition for different modes 
and also impose that the time components vanish at the horizon, 



(47) 



where we have used the linearity of eqs. (45) to set the scale of at the horizon to 1. The 
coefficients in the above expansions can be fully determined by plugging the expansion into eqs. (45) 



4 


= (1 


-u) 


-i£>/4 


) + 4 {1 \± 


-u)+af\l 


-uf + 4 i3) (i 


-uf + -..' 


a] 


= (1 


-«)■ 


-iui/A 




+ af\l-u) 




3 + ..." 


5 


4 


= (1 


-«)- 


-iui/A 






2 + af\l-u 


3 + ..." 


; 



and counting powers of (1 — u). Then, eqs. (47) can provide initial conditions for these second 
order differential equations (fi"5|). In fact, we use these power solutions to do numerical integration 
from the horizon to the conformal boundary as mentioned in the last section. 
At the conformal boundary, the general solutions to eqs. (Ha) are of the form, 



,1(0) 



,i(i) 



A l t " + A$ K '\ogu + --- 



£ = A 



2(0) 



A 2(1) logn + 



(48) 



We can also expand the constraint equations near the conformal boundary, but it is of no use for 
later calculations. 

As argued in 0], for the conductivity to be a gauge invariant quantity, we need to construct a 
new mode from a|,a 2 ,a^ which should be invariant under the gauge transformation that respects 



7 These constraints come form the radial gauge we have chosen for the gauge field fluctuation, i.e., n° = when 
deriving the equations of motion. 
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our gauge choice. We here do not go into the details of the construction of this field but write 
down the mode directly, 



iua? + 6a} 



al = a% + i> J T \ 49 



Plugging the boundary behavior as in eqs. (48) into the newly defined mode and expanding it near 
the conformal boundary, 

al = Af ) +Af ) u 2 + ..., (50) 

with 

Af' = Af\ Ai m = ^> + Af\ (51) 

A 4 

Then, the conductivity along the x direction is straightforwardly defined as 

1 i3 (1) 

Vxx(v) = ^T^V- (52) 

v ; iu 7 3 (°) V ; 

Notice that the above formula is very different from equation (4.19) of Q due to the appearance 
of the logarithmic terms in the time component of the SU(2) gauge field (for both the background 
and the fluctuations). The definition for the spectral function is similar to the s-wave case and we 
list them in what follows for completeness, 

n = -2%y (53) 

along the x-direction; 

U = - 2 ^7nT ( 54 ) 

along the y-direction. 

In FIG. [ll] we plot numerical results for the conductivity along the x direction. We can see 
that the results are very similar to those of the y direction stated in FIG. pi]. Explicitly, there 
is a DC infinity in the conductivity along the x direction. What is strikingly different from the 
AdSi p-wave superconductor is that there is no pole for the imaginary part of the conductivity at 
nonzero frequency. Mathematically, this is due to the logarithmic term appearing in the boundary 
behavior for the time component fluctuations of the gauge field. We will take this as the effect 
of the nontrivial dynamical exponent z / 1, The real (imaginary) part of the conductivity ap- 
proaches 1 (0) more slowly than the y direction case, which is consistent with the spectral function 
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FIG. 11. The real (imaginary) part of the AC conductivity along x direction of the p-wave superconductor 
at T/T c = 1.0,0.588587,0.171829,0.079889,0.0691493 from top to bottom (bottom to top). 




to 



T 

FIG. 12. The spectral functions for the modes a 3 y (red) and (blue) at temperature corresponding to the 
conductivity pictures. 



profiles plotted in FIG. 12. This can be explained by the anisotropic characteristic of the p-wave 
superconductor. 

The qualitative behavior is similar to the corresponding s-wave case: at high frequency the real 
part of the conductivity approaches 1 while the imaginary part goes to 0; a peak also appears at 
zero frequency, indicating a DC superconductivity; the gap forms as the temperature decreases and 
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it gets increasingly deep until the conductivity is exponentially small. The appearance of the delta 
peak at zero frequency can also be well understood by the Kramers-Kronig relation ( j35|) . Although 
the conductivity a yy (uj) behaves qualitatively similar to ct xx (uj), the anisotropic characteristic of 
the p-wave superconductor can also be seen from the two conductivities: the real part of the (t xx (oj) 
grows much slowly than that of a yy (u). These features are common with the AdS superconductors. 
One main difference is that there is no pole at nonzero frequency for the imaginary part of a xx (u) 
as mentioned before. Another main difference from the AdS black hole superconductor is that the 
imaginary part of the conductivity never goes below zero and approaches zero quite slowly. 



V. SUMMARY 



In this work we explored properties of holographic superconductors with nontrivial dynamical 

exponent by putting the Abelian-Higgs model (s-wave) or SU(2) gauge field (p-wave) into the 

Lifshitz black hole geometry constructed in ||. We found that the order parameters (O) all have 

mean field behavior (T — T c ) 1 / 2 near the critical temperature T c , which is qualitatively consistent 

with the AdS superconductors as well as BCS theory. One difference between our results and 

previous investigations on condensates is that the condensates do not approach some fixed values 

in the zero temperature limitQ. We then plot the free energy difference between the normal and 

superconducting phases, which can be taken as evidence of the occurrence of the superconducting 

phase transition. We also numerically compute the AC conductivities and they nearly behave in 

the same way: a peak appears at zero frequency, indicating a DC superconductivity; the gap forms 

as the temperature decreases and it gets increasingly deep until the real part of the conductivity 

gets exponentially small. The anisotropic characteristic of the p-wave superconductor can be seen 

from the difference between a xx (u) and a yy {uj): their growth rates as increasing the frequency are 

very different. But we do not see a pole at nonzero frequency for the imaginary part of a xx (uj), 

which does exist in the p-wave superconductor when taking the AdS^ black hole geometry. Another 

feature is that all the real parts of the conductivities approach 1 (but never exceed this value) in 

contrast to those of the AdS case. With respect to this, the imaginary parts of the conductivities 

approach zero in the high frequency limit but never go below zero. We attribute these differences 

from the AdS case to the effect of nontrivial dynamical exponent. More specifically, the black 

hole geometry considered in this work is anisotropic between space and time, very different from 

8 We would like to state that this feature is not due to the numerical approach used here. Actually, we already used 

our numerical method to carry out similar computations in other black hole backgrounds pi] ], of asymptotic AdS 

type and found that the condensed operator approach constant value. Therefore, we conclude that the feature 

found here is due to the nontrival dynamical exponent z. 

22 



the Schwarzschild-AdS black hole, which results in different asymptotic behaviors of temporal and 
spatial components of gauge fields than previous conclusions in Schwarzschild-AdS black hole. 
These common features also imply that general gauge/gravity duality is a useful and powerful tool 
in producing some universal properties of strongly coupled system in condensed matter physics. 



ACKNOWLEDGMENTS 

The author would like to thank Johanna Erdmenger, Jonathan Shock, Xin Gao and Da- Wei 
Pang, Xu Zhang for useful discussions. This work was supported by the MPS-CAS Doctoral 
Training program. 



[1] J. M. Maldacena, |Adv.Thcor.Math.PhysJ 2, 231 (1998), |arXiv:hcp-th/9711200 [hcp-th] 



[2] J. Erdmenger, N. Evans, I. Kirsch, and E. Threlfall, |Eur.Phys.J.| A35, 81 (2008), |arXiv:071 1.4467 [hep- 



_th]|S. S. Gubser and A. Karen, |Ann.Rcv.Nucl.Part.SciJ 59, 145 (2009), |arXiv:0901.0935 [hcp-th]| 5. A. 
Hartnoll, |Class.Quaiit.GravJ 26, 224002 (2009), |arXiv:0903.3246 [hep-thl c. P. Herzog, [j.Phys.A| A42, 
343001 (2009), |arXiv:0904.1975 [hep-thf U. Gursoy, |Mod.Phys.LcttJ A23, 3349 (2009), |arXiv:0904.2750 



[hep-th| J. McGreevy, |Adv.High Energy Phys.| 2010, 723105 (2010), |arXiv:0909.0518 [hcp-th] 



[3] S. A. Hartnoll, C. P. Herzog, and G. T. Horowitz, |Phys.Rev.Lett.| 101, 031601 (2008), |arXiv:0803. 3295 
[hcp-th|jHEP| 0812, 015 (2008), |arXiv:0810.1563 [hcp-th 



[4] S. S. Gubser and S. S. Pufu, [JHEP| 0811, 033 (2008), |arXiv:0805.2960 [hcp-th 



[5] S. Kachru, X. Liu, and M. Mulligan, |Phys.Rcv.| D78, 106005 (2008), |arXiv:0808.1725 [hcp-th] 



[6] G. Bertoldi, B. A. Burrington, and A. Peet, |Phys.Rcv.| D80, 126003 (2009), |arXiv:0905.3183 [hcp- 
th] | M. Taylor(2008), brXiv:0812.0530 [hcp-th] 



[7] P. Kovtun and D. Nickel, [Phys.RevLettJ, 102, 011602 (2009), brXiv:0809.2020 [hep-th] D.-W. Pang, 
|JHEP| 1001, 120 (2010), fc.rXiv:0912.2403 [hcp-th] E. Brynjolfsson, U. Danielsson, L. Thorlacius, and 
T. Zingg, 1008, 027 (2010), |arXiv:1003.5361 [hep-th| C. Hoyos-Badajoz, A. O'Bannon, and J. M. 
Wu, 1009, 086 (2010), |arXiv: 1007.0590 [hep-tl^| T. Nishioka and H. Tanaka, 1102, 023 (2011), 



arXiv: 1010.6075 [hcp-th| S. A. Hartnoll, D. M. Hofman, and D. Vegh, 1108, 096 (2011), |arXiv:1105.3197 



[hep-thp . P. Lcmos and D.-W. Pang, 1106, 122 (2011), |arXiv:1106.2291 [hcp-th]| L. Huijse, S. Sachdev, 



and B. Swingle, |Phys.Rcv.| B85, 035121 (2012), |arXiv:1112.0573 [cond-mat.str-cl]| Y. S. Myung, 
|Eur.Phys.J.| C72, 2116 (2012), bjXiv:1203.1367 [hcp-th] 



S. A. Hartnoll, J. Polchinski, E. Silverstein, and D. Tong, |JHEP| 1004, 120 (2010), |arXiv:0912.Tu61 



[hep-th] 



T. Azeyanagi, W. Li, and T. Takayanagi, |JHEP| 0906, 084 (2009), |arXiv:0905.0688 [hcp-th 



23 



[10] L. Q. Fang, X.-H. Ge, and X.-M. Kuang(2012), |arXiv:1201.3832 [hep-th]| M. Alishahiha, M. Moham- 
madi Mozaffar, and A. Mollabashi(2012), |arXiv:1201.1764 [hcp-th] 



[11] S.-J. Sin, S.-S. Xu, and Y. Zhou, Int.J.Mod.Phys. A26, 4617 (2011), arXiv:0909.4857 [hep-th 



[12] E. Brynjolfsson, U. Danielsson, L. Thorlacius, and T. Zingg, J.Phys.A A43, 065401 (2010), 



arXiv:0908.2611 [hep-th 



[13] R.-G. Cai and H.-Q. Zhang, Phys.Rev. D81, 066003 (2010), arXiv:091 1.4867 [hep-th 



[14] U. H. Danielsson and L. Thorlacius, JHEP 0903, 070 (2009), arXiv:0812.5088 [hcp-th] R. B. Mann, 



0906, 075 (2009), |arXiv:0905.1136 [hcp-th][ D.-W. Pang(2009), |arXiv:0905.2678 [hcp-th| D.-W. Pang, 
0910, 031 (2009), pjXiv:0908.1272 [hcp-thjl K. Balasubramanian and J. McGreevy, |Phys.Rev.| D80, 
104039 (2009), |arXiv:0909.0263 [hcp-th"l E. Ayon-Bcato, A. Garbarz, G. Giribct, and M. Hassainc, 
D80, 104029 (2009), ^rXiv:0909.1347 [hcp-th]| R.-G. Cai, Y. Liu, and Y.-W. Sun, |JHEP| 0910, 080 
(2009), |arXiv:0909.2807[hep-th1| Y. S. Myung, Y.-W. Kim, and Y.-J. Park, |Eur.Phys.j] C70, 335 (2010), 



|arXiv:0910.4428 [hep-thj D.-W. Pang, |JHEP| 1001, 116 (2010), |arXiv:0911.2777 [hcp-th]| K. Goldstein, 
S. Kachru, S. Prakash, and S. P. Trivedi, 1008, 078 (2010), |arXiv:091 1.3586 [hcp-th| E. Ayon-Beato, 
A. Garbarz, G. Giribet, and M. Hassaine, 1004, 030 (2010), |arXiv:1001.2361 [hcp-thf j. Blaaback, U. H. 
Danielsson, and T. Van Rict, 1002, 095 (2010), |uXiv:1001.4945 [hcp-th| Y. S. Myung, |Phys.Lctt.| B690, 
534 (2010), |a,rXiv:1002.4448 [hep-th] 



[15] C.-M. Chen and D.-W. Pang, |JHEF| 1006, 093 (2010), |arXiv: 1003.5064 [hep-thj| K. Balasubrama- 
nian and K. Narayan, 1008, 014 (2010), |arXiv:1005.3291 [hcp-th]|K. Goldstein, N. Iizuka, S. Kachru, 



S. Prakash, S. P. Trivedi, et al, 1010, 027 (2010), |arXiv: 1007.2490 [hep-th]| A. Donos and J. P. Gauntlett, 
1012, 002 (2010), ^.rXiv: 1008.2062 [hcp-thj| R. Gregory, S. L. Paramcswaran, G. Tasinato, and I. Zavala, 



1012, 047 (2010), arXiv: 1009.3445 [hcp-th] M. Dchghani, R. Mann, and R. Pourhasan, Phys.Rev 



D84, 046002 (2011), |arXiv:1102.0578 [hep-t,h"]] D. Cassani and A. F. Faedo, [JHEP| 1105, 013 (2011), 



arXiv:1102.5344 [hcp-th| I. Amado and A. F. Faedo, 1107, 004 (2011), |arXiv:1105.4862 [hcp-th]| j. Tarrio 



and S. Vandorcn, 1109, 017 (2011), |uXiv:1105.6335 [hep-th]^ . F. Ross, |Class.Quant.GravJ 28, 215019 
(2011), |arXiv:1107.4451 [hep-trif x. Dong, S. Harrison, S. Kachru, G. Torroba, and H. Wang, |JHEP 
1206, 041 (2012), |arXiv:1201.1905 [hcp-thf Y. Liu, 1206, 024 (2012), |arXiv:1202.1748 [hcp-th| L. Bar- 



clay, R. Gregory, S. Parameswaran, G. Tasinato, and I. Zavala, 1205, 122 (2012), arXiv:1203.0576 



[hep-th] 



[16] S. S. Gubser, Phys.Rev. D78, 065034 (2008), arXiv:0801.2977 [hcp-th 



[17] Y. Bu(2012),|arXiv:1205.1614 [hep-th] 



[18] I. R. Klebanov and E. Witten, |Nucl.Phys.| B556, 89 (1999), |arXiv:hcp-th/9905104 [hcp-th] 



[19] D. T. Son and A. O. Starincts, JHEP 0209, 042 (2002), arXiv:hcp-th/0205051 [hcp-th] 



[20] M. Amnion, J. Erdmenger, M. Kaminski, and P. Kerner, |Phys.Lettj B680, 516 (2009), |arXiv:0810. 2316 



[hcp-th|JHEP| 0910, 067 (2009) , |arXiv:0903. 1864 [hcp-th][ P. Basu, J. He, A. Mukherjee, andH.-H. Shieh, 



0911, 070 (2009), ^rXiv:0810.3970 [hcp-th| K. Peeters, J. Powell, and M. Zamaklar, 0909, 101 (2009), 



arXiv:0907.1508 [hep-th 



24 



[21] Y.-Y. Bu, |Nucl.PhysJ B864, 806 (2012) 



25 



